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Single-particle diagonal and off-diagonal Green's functions of a 2-leg t-J ladder at 1/8-doping 
are investigated by Exact Diagonalisations techniques. A numerically tractable expression for the 
superconducting gap is proposed and the frequency dependence of the real and imaginary parts of 
the gap are determined. The role of the low-energy gapped spin modes whose energies are computed 
by a (one-step) Contractor Renormalization procedure are discussed. 
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In the BCS theory of superconductivity [l| the fre- 
quency and momentum dependence of the superconduct- 
ing gap provide information about the frequency and 
momentum dependence of the pairing interaction. In 
one spatial dimension quantum phase fluctuations de- 
stroy long-ranged SC order although the spectral gap is 
expected to survive 2]. In a doped two-leg spin lad- 
der Q the doped holes can form mobile singlet pairs Q 
leading to dominant algebraic SC correlations |5| and a 
robust spin gap. In this Letter, motivated by previous 
work by two of the authors 0, we introduce a numer- 
ically tractable expression for the superconducting gap 
function which, as we argue, contains a non-trivial mo- 
mentum and frequency dependence apart from an overall 
(vanishing in the thermodynamic limit) prefactor which 
accounts for ID quantum phase fluctuations. 

We shall consider here a generic 2-leg t-J ladder, 

H = Jleg <Sj,o • Sj+l,a + -/rung Sj t i • Sjfi (1) 
i,a i 

+ *leg 4,0^+1.0 + *rung ^ 4,a C M+l + H C - ' 
i,a i,a 

where Ci^ a are projected hole operators and a (=1,2) 
labels the two legs of the ladder. Isotropic couplings, 



t and Jic K = J n 



J, will be of interest 



here. A value like J = 0.4 (t is set to 1) and a doping 
of 6 = 1/8 are typical of superconducting ladder materi- 
als such as Sri4_ a; Ca a ;Cu2404i and these parameters will 
be assumed hereafter. Single-particle diagonal and off- 
diagonal spectral functions are computed by Exact Diag- 
onalisations of 2 x L periodic ladders of size L — 12. The 
expression for the superconducting pairing function in- 
troduced previously in the context of the two-dimensional 
tJ model 6] is computed numerically from the knowl- 
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FIG. 1: Low-energy magnetic excitations at 1/8-doping and 
J — 0.4 for q y — (circles) and q y = ir (diamonds) mo- 
menta, (a) Triplet collective mode obtained by CORE on 
2 x 16 (shaded symbols) and 2 x 24 (closed symbols) ladders. 
The particle-hole continuum (open symbols) constructed from 
the data of Figs. |5|&|S|is also shown; (b) S = 3/2 excitations 
(measured from the chemical potential) calculated by ED of 
the 2 x 12 cluster. The onsets of the 3-quasi-particle continua 
are shown by arrows. 



edge of the diagonal and off-diagonal Green's functions 
of the doped 2-leg ladder. The role of the low-energy 
collective triplet modes whose energies are obtained by a 
Contrator-Renormalisation 0, (CORE) calculation is 
investigated. 

As established by numerical or Bosonisation tech- 
niques, in the parameter range considered here, the 2- 
leg doped spin ladder exhibits dominant superconducting 
fluctuations 0] and its low-energy physics is governed by 
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FIG. 2: Spectral function vs u) (t = 1) on a 2 x 12 ladder at FIG. 3: Momentum dispersion of the low-energy peaks seen 

1/8-doping and J = 0.4 for two transverse momenta and in Fig.|2] Error bars give the widths of the peaks (if any) and 

7r in panels (a) and (b) respectively. Data for different chain the dotted lines correspond to the lowest excitation energies, 
momenta are shifted w.r.t each other. 



two (weakly gapped) collective spin modes (magnons), 
a gapped (q y = it) charge mode and a (q y — 0) zero- 
energy collective charge mode 0, 0, Q| characteristic of 
a C1S0 phase of the Luther-Emery (LE) liquid universal- 
ity class ^3 . Prior to the investigation of spectral prop- 
erties we have computed the lowest magnon excitations 
at hole density 1/8 on ladders with up to size 2 x 24 with 
Nit = 6 holes using an effective CORE hamiltonian [j| as 
shown on Fig. ^a) . Our extrapolation gives a spin gap 
(for q y = 7r) AJ^g ~ 0.11 significantly smaller than the 
spin gap ~ J/2 = 0.2 of the undoped spin ladder. On the 
other hand, the lowest q y — triplet (not shown) occurs 
at A° lag ~ 0.17, close to the onset of the particle-hole 
continuum (discussed below). 

We now turn to the investigation of the single-particle 
excitations. For a finite size system, it is convenient to 
define Green's functions so that the sets of electron-like 
(for oj < 0) and hole-like (for uj > 0) poles coincide. The 
time-ordered diagonal Green's function is then defined 
as, 



G(q,w) = (N\c\ 



uj — ie + H — £7jv-i 

-— )— cl a \N-2), (2) 

uj + ie — H + Ln-i 

and the retarded one, G re t(q, oj), by the substitution 
u + ie — > oj — ie in the second term, e is a small imaginary 
part set to 0.05 hereafter. Note that the "symmetriza- 
tion" between oj and — oj implies using two ground states 
(GS) with particle numbers N and N — 2 surrounding 
the value N — I = N s — Nh corresponding to the actual 
hole density ie here Nh = 3 holes on N s = 2L = 24 sites. 
For convenience, the energy reference -Ejv— i is defined as 
the average between the GS energies of N} and |iV — 2). 
G(q, lu) and G rc t (q, w) can be calculated using a standard 



continued-fraction method. 

The data for the single-particle spectral function 
A(q, oj) = — ilmG re t(q, shown in Fig. |2]can be fairly 
well described by (i) BCS-like q y — (bonding) and 
q y = n (antibonding) quasi-particle bands with the chem- 
ical potential reference energy u> = and (ii) a broad 
incoherent background extending further away towards 
negative energies in agreement with a calculation on a 
smaller 2x8 cluster at the same hole density 13] . Note 
that the peaks of the low-energy band-like features nar- 
row near the chemical potential. Note also that the ap- 
proximate Fermi momenta (feF,i>0) which lies between 
(57r/6,0) and (2n/3, 0) and (fci?2,7r) ~ (7r/3, it) are in 
rough agreement with Luttinger's theorem | l4| (which 
gives kpi+kp2 = 7tt/8). The low-energy peaks of A(q, uj) 
plotted vs momentum in Fig. [3] exhibit BCS-like disper- 
sions ±^je 2 ^ + A 2 (g) with the magnitudes of the gaps 

A(fc F i,0) = A ~ 0.115 and A(fc F2 ,7r) = A,, ~ 0.090. 
Here, the solid points denote the maximum spectral 
weight and the open circles the BCS-like shadow band. 
Note that the particle-hole continuum can been obtained 
by considering all combinations of any two of the lowest 
single-particle excitations. Figs, ^a) and^b) give fur- 
ther evidence that bound triplet-hole pair excitations sits 
below this continuum j9j. 

The superconducting Gorkov's off-diagonal one- 
electron time-ordered Green's function Q] defined by 



F(q,i) = »(Tc_ q ,_ (r (t/2)c q , <r (-*/2)> , 



(3) 



can be computed in a finite system by taking the expecta- 
tion value between the two GS |iV) and \N — 2) differing 
by two particles, hence reflecting superconducting fluc- 
tuations. Note that these states are both spin singlets so 
that the (Fourier transformed in time) Green-function, 
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FIG. 4: Spectral weight of the superconducting Green's func- 
tion. Same conventions as in Fig. |5] Note the change of sign 
between transverse momenta and n. 



.F(q, uj), is even in frequency and can be expressed as, 

F(q, w) = # q (w + ie) + F q (-w + te) , (4) 



where 



1 



F q (z) = (TV - 2|c_ q ,_ ff — , c q , CT |iV>, (5) 

Z — 11 + hiN-l 

is defined for all complex z (with Imz ^ 0). Follow- 
ing Ohta, et. al [l5j] who extended the continued-fraction 
method to deal with off-diagonal Green's functions, we 
have computed F(q, uj) on a 2x 12 ladder at a hole density 
of 1/8 and for all momenta. Sharp low energies features 
are found near momenta (/cf,i,0) and (fcF,2)7r) as shown 
in Fig. 2| As the diagonal Green's function, the off- 
diagonal one is also gapped. The equal-time pair ampli- 
tude (c_ q) _ CT Cq >CT ) corresponding to the integrated weight 
I^oo ^Xl) w ) ^ ^ s plotted in Fig. [3] showing maxima 
in the vicinity of the (estimated) Fermi momenta. Note 
also in Figs. 0] and [S] that the pair amplitudes have oppo- 
site signs for the two transverse momenta and 7r, rem- 
iniscent of a two-dimensional (2D) d-wave orbital sym- 
metry. 

A long-ranged superconducting GS is characterized by 
a frequency-dependent gap function Asc (q, uj) directly 
proportional to F(q, u>). A definition such as 



A S c(q, uj) 



2uF{q,uj) 



(6) 



G(q,w) - G(q,-u) ' 

is consistent with Nambu-Eliashberg theory as shown 
by two of the authors |6j. From a numerical calculation 
of G(q, uj) and F(q,u>), Asc(q, uj) has been obtained in 
a 2D t-J model extending previous work by Ohta 
et. al ^3 who fit the spectral weight ImF(q, w)/7r to a 
d x 2_ y 2 BCS-Bogoliubov quasi-particle form [lj. In con- 
trast to a true SC state, a doped 2-leg ladder (for param- 
eters used here) exhibits long distance power law pair 



field correlations which decay as x~ 1 ^ Kp . Here K p is 
the Luttinger liquid parameter associated with the mass- 
less charge mode. This implies 0] that for a ladder of 
length L, the off-diagonal Green's function F(q,u>) de- 
cays as (£/ 'L) 1 / 2K f> . Here the coherence length £ is pro- 
portional to the inverse of the spin gap. Thus, we expect 
that A(q,w) given by Eq. © will vary as (^/L) 1 / 2K p. 
Using a CORE calculation supplemented by conformal 
invariance identities ||, we obtain, for S = 1/8 and 
J = 0.4, K„ ~ 0.65 in agreement with previous ED eval- 
uations [13] and DMRG data [13 • We then expect the 
SC gap function Asc(q, uj) to vanish in the thermody- 
namic limit due to Cooper pair phase fluctuations. In- 
deed, using e.g. a standard low-energy long-wavelength 

LE field theory, it can be shown that the SC Green's 

i 

function decays with system size like (1/L) 2K p due to SC 
phase fluctuations However, apart from this pref ac- 
tor, Asc(q ; uj) calculated on a finite system can provide 
informations on the dynamics of pairing at intermediate 
distances. 

In our case L/£ is only of order 2 to 3 so that the scaling 
factor (£/L) 1 / 2Kp is of order one and we will simply nor- 
malize 18] As C (q,u) so that A SG (k F1 , 0, uj = 0) = 0.12/ 
and Asc(/cp2, i", uj = 0) = —0.09/. Using this normaliza- 
tion, we have plotted the real and imaginary parts of 
Asc(qjw) for the bonding q = (27r/3,0) and antibond- 
ing q = (tt/3, 7t) Fermi momenta in Fig. As expected 
for a d-wave-like gap, the sign of the gap changes when 
one goes from (kpi,0) to (&F2,7r). Otherwise, the fre- 
quency dependence of both the real and imaginary parts 
of the gap is quite similar. The imaginary part of the 
gap appears to onset at values of ui ~ Ao + AJj ag where 
Ao ~ 0.12 is the superconducting gap and A^ ag ~ 0.11 
is the magnon gap. The imaginary part of the gap then 
increases until one passes through the particle-hole spec- 
trum and then decreases at yet higher energies. The 
real part of the gap increases slightly and then when uj in- 
creases beyond the electron-hole spectrum, the real part 
of the gap drops. 

Thus, we believe that this approach to calculating 
Asc(qjw) allows one to probe the internal structure of 
a pair. Clearly, additional calculations, particularly for 
the 2-leg Hubbard ladder will be of interest in providing 
further insight into the relationship between frequency 
dependence of the gap and the dynamics of the pairing 
interaction. 
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FIG. 5: Superconducting amplitude as a function of the chain 
momentum q x for the bonding q y = (solid circles) and an- 
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FIG. 6: Real (a) and imaginary (b) parts of the superconduct- 
ing gap vs ui, normalization as discussed in the text, obtained 
by ED results for a 2 x 12 t-J ladder at the two bonding 
and antibonding Fermi momenta (27r/3, 0) and (7r/3, n). An 
imaginary damping e = 0.05 is used and the irrelevant offset 
F(0, 0) oc ie has been subtracted. 
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